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Summary 

Epistasis is an important feature of the genetic architecture of quantitative traits, but the dynamics of epistatic 
interactions in natural populations and the relationship between epistasis and pleiotropy remain poorly 
understood. Here, we studied the effects of epistatic modifiers that segregate in a wild-derived Drosophila 
melanogaster population on the mutational effects of P-elenient insertions in Semaphorin-SC {Sema-Sc) and 
Calreticulin (Crc), pleiotropic genes that affect olfactory behaviour and startle behaviour and, in the case of Crc, 
sleep phenotypes. We introduced Canton-S B (CSB) third chromosomes with or without a /"-element insertion at 
the Crc or Sema-5c locus in multiple wild-derived inbred lines of the Drosophila melanogaster Genetic Reference 
Panel (DGRP) and assessed the effects of epistasis on the olfactory response to benzaldehyde and, for Crc, also 
on sleep. In each case, we found substantial epistasis and significant variation in the magnitude of epistasis. The 
predominant direction of epistatic effects was to suppress the mutant phenotype. These observations support a 
previous study on startle behaviour using the same D. melanogaster chromosome substitution lines, which 
concluded that suppressing epistasis may buffer the effects of new mutations. However, epistatic effects are not 
correlated among the different phenotypes. Thus, suppressing epistasis appears to be a pervasive general feature 
of natural populations to protect against the effects of new mutations, but different epistatic interactions 
modulate different phenotypes affected by mutations at the same pleiotropic gene. 



1. Introduction 

Epistasis is an integral feature of the genetic 
architecture of quantitative traits (Anholt & Mackay, 
2004; Flint & Mackay, 2009; Mackay et ai, 2009). 
Epistasis occurs when the effect of variation at one 
locus is suppressed or enhanced by the genotype at 
another locus. Epistatic interactions can bias esti- 
mates of the effects of quantitative trait loci (QTLs) in 
mapping populations when present but not accounted 
for (Carlborg et al., 2006) ; enable inferences of genetic 
networks affecting complex traits (Phillips, 2008) ; and 
affect predictions of long-term response to artificial 
and natural selection (Carlborg et al., 2006; Phillips, 
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2008). Epistasis is difBcult to detect in classical quan- 
titative genetic analyses based on resemblance be- 
tween relatives in outbred populations (Falconer & 
Mackay, 1996), and epistatic interactions contribute 
largely additive genetic variation in outbred popu- 
lations when the contributing alleles are rare (Hill 
et al., 2008). However, epistatic interactions are com- 
mon in experiments designed to examine their effects 
on trait means in QTL mapping populations. For 
example, in Drosophila epistatic interactions have been 
reported between QTLs affecting bristle number (Long 
et al., 1995; Gurganus et al., 1999; Dilda & Mackay, 
2002), wing morphology (Weber et al., 1999), lifespan 
(Leips & Mackay, 2000, 2002) and startle-induced 
locomotor behaviour (Jordan et al., 2006). In mice, 
epistasis has been reported between QTLs affecting 
growth, body weight and morphometry (Brockmann 
et al., 2000; Cheverud et al., 2001; Workman et al., 
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2002; Klingenberg et al., 2004; Yi et al., 2006). 
Epistasis is also a prominent feature of the genetic 
architecture of growth rate in Arabidopsis (Kroymann 
& Mitchell-Olds, 2005), chickens (Carlborg et al., 
2006) and yeast (Steinmetz et al., 2002; Sinha et al., 
2008). 

Although epistatic interactions have been detected 
in QTL mapping experiments, it is easier to study 
epistasis in crosses among Hues with reduced genetic 
heterogeneity in largely homozygous genetic back- 
grounds (Eshed & Zamir, 1996; Clark & Wang, 1997; 
Sambandan et al., 2006). Drosophila melanogaster is 
an excellent model system to study epistasis affecting 
quantitative traits due to the ease of constructing 
chromosome substitution and introgression lines, and 
generating mutations in a common homozygous 
genotype. Epistasis has been documented for ag- 
gressive behaviour by constructing chromosome 
substitution Unes in which small segments of one 
genotype were introgressed into a different genetic 
background (Edwards & Mackay, 2009). Epistasis for 
aggression was also evident from behavioural and 
whole genome transcriptional analyses of an ensemble 
of co-isogenic hyper-aggressive _P-element insertion 
hues (Zwarts et al., 2011). Epistasis for metabolic 
activity was revealed by constructing aU possible 
two-locus genotypes for several pairs of /"-element 
insertion mutations (Clark & Wang, 1997). Diallel 
cross analysis of co-isogenic P-element insertion 
hues enabled identification of epistatic networks of 
genes affecting negative geotaxis (Van Swinderen & 
Greenspan, 2005), olfactory avoidance behaviour 
(Fedorowicz et al., 1998; Sambandan et al., 2006), 
aggression (Zwarts et al., 201 1) and startle behaviour 
in Drosophila (Yamamoto et al., 2009). 

Previously, Yamamoto et al. (2009) created pairs 
of chromosome substitution lines in which isogenic 
CantonS B (CSB) chromosomes with /"-element in- 
sertions in genes affecting startle behaviour and their 
/"-element free co-isogenic control chromosomes were 
substituted into different homozygous wild-derived 
D. melanogaster genotypes. This design enables the 
quantification of the extent to which naturally segre- 
gating variants modify the effects of single mutations, 
as well as the magnitude of variation of epistasis 
among the different lines. This study reported wide- 
spread suppressing epistasis of naturally segregating 
modifiers on startle behaviour. Since the magnitude of 
suppressing epistasis was proportional to the magni- 
tude of the mutational effect of the /"-element inser- 
tion on startle behaviour, it was concluded that 
suppressing epistasis buffers the effects of new muta- 
tions in natural populations. 

P-element insertions at genes previously implicated 
in startle behaviour, Seniaphorin-SC {Senw-Sc) and 
Calreticulin (Crc) also affect olfactory behaviour 
(Sambandan et al., 2006) and in the case of Crc, sleep 



phenotypes (Harbison & Sehgal, 2008). The objective 
of the present study was to ask whether suppressing 
epistasis by naturally segregating modifiers on behav- 
ioural traits is a general principle or unique to the 
startle response, and, moreover, to assess whether the 
effects of the same P-element insertion on different 
phenotypes is modulated by the same or different 
epistatic modifiers. 

2. Materials and methods 

(i) Drosophila stocks 

/"-element insertion fines for Crc and Sema-Sc, which 
were generated as part of the Berkeley Drosophila 
Gene Disruption Project (Bellen et al., 2004), contain 
single P[GT1] insertions generated in the isogenic 
u'^^^*, CSB background. Crc and Sema-5c have pleio- 
tropic effects on olfactory avoidance of benzaldehyde 
(Sambandan et al., 2006; Rollmann et al., 2007), 
bristle number (Norga et al., 2003) and startle re- 
sponse (Yamamoto et al., 2008, 2009), and Crc also 
has pleiotropic effects on sleep traits (Harbison & 
Sehgal, 2008). Both Crc and Sema-5c are located on 
chromosome 3 (Ci). The construction of chromo- 
some substitution lines carrying either a CSB C3 or 
the Crc and Sema-5c P[GT1] mutations on the same 
CSB Ci in inbred lines of the Drosophila melanogaster 
Genetic Reference Panel (DGRP ; Mackay etal., 2012), 
derived from a Raleigh (North Carohna) population 
of wild D. melanogaster, has been reported previously 
(Yamamoto et al., 2009). Thirteen chromosome sub- 
stitution fines with /"-element insertions at Sema-5c 
and 14 chromosome substitution fines with /"-element 
insertions at Crc and the corresponding controls were 
used in this study (Fig. 1). AU flies were reared in large 
mass cultures on cornmeal/molasses/agar medium at 
25 "C and a 12 h light/12 h dark cycle (fights on at 
6 : 00 am; lights off' at 6 : 00 pm). 

(ii) Behavioural assays 

We measured olfactory behaviour for CJ substitution 
lines with Crc and Sema-5c mutations, and the cor- 
responding CSB Ci substitution fines contempo- 
raneously using a modification of the ' dipstick ' assay 
(Anholt et al., 1996), as described previously (Swarup 
et al., 2011). We measured olfactory behaviour 
of single-sex groups of 50 flies/replicate and three 
replicates/sex for each fine. Assays were conducted 
between 2 : 00 and 4 : 00 pm using 0-3 % (v/v) benzal- 
dehyde (Sigma-Aldrich, St. Louis, MO). Replicate 
measurements on individual fines were collected over 
multiple days to account for environmental variation. 
Flies between 4 and 7 days old were coUected a day 
prior to the assay and food deprived for 2 h in a 50 ml 
conical tube containing a cotton wool swab tip 
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Fig. 1 . Generation of co-isogenic CSB C3 substitution lines in inbred DGRP genetic backgrounds. Ttie left side of 
the diagram illustrates the three major D. melanogaster chromosomes in co-isogenic CSB lines, with arrows indicating 
the locations of P-element insertions in Sema-Sc and Crc. The right side of the diagram illustrates the introduction of 
CSB C3 with or without P-element insertions into different DGRP Unes, indicated with different colours (Yamamoto et al., 
2009). 



(referred to as 'odour tube'). The measurement is 
initiated by depositing 01 ml of odorant solution on 
the cotton wool swab tip in the odour tube. The odour 
tube is then connected to a collection tube and flies are 
given 2 min to partition between the tubes. At the end 
of the assay, a response index (RI) is calculated as 
follows: Rl = number of flies in the collection tube/ 
total number of flies. An Rl of 1 indicates the highest 
avoidance response to benzaldehyde, while 0 indicates 
indifference (or attraction) to the odorant. 

We measured sleep and waking activity of the Crc 
chromosome substitution lines and their respective 
controls by recording locomotion of virgin male 
and female flies for seven continuous days using the 
Drosopliila Activity Monitoring (DAM) System 
(Trikinetics, Waltham, MA). Each fly was housed 
separately in an activity monitor tube. The DAM sys- 
tem uses an infrared beam to detect movement in the 
monitor tube; the movement is recorded as activity 
counts in 1-min intervals. We eliminated flies that died 
after 7 days of recording from the analysis. We used a 
custom C++ program to compute day and night sleep 
duration in minutes, and waking activity as counts per 
waking minute. Sleep is defined as inactivity lasting 
5 min or longer (Hendricks et al., 2000; Shaw et al., 
2000; Huber et al., 2004, Ho & Sehgal, 2005). 

(iii) Mutational effects and epistatic interactions 

We estimated the effects (2a) of each mutation on ol- 
factory behaviour and sleep phenotypes in the CSB 
background as the deviation of the mean phenotypic 
value of the homozygous mutant from that of the 
CSB control (Falconer & Mackay, 1996). We used 



Student's t tests to assess the significance of the dif- 
ference in phenotypic values between mutant and 
control. 

We estimated the epistatic interaction for each 
DGRP line as the difference between the expected and 
observed phenotypic values. There are two chromo- 
some substitution fines for each DGRP line, one with 
a mutant C3 and the other with a wild-type CJ. The 
observed phenotypic value of each DGRP fine is tfie 
mean of the fine with the mutant Ci. The expected 
phenotypic value of each fine is the difference between 
the mean of the line with the wild-type Ci and the 
estimate of 2a for the appropriate mutation obtained 
in the pure CSB background. We assessed the sig- 
nificance of epistatic interactions in each DGRP line 
background by performing tfiree-way fixed effect 
analyses of variance (ANOVAs) using the model : Y= 
pi+G + L + S+GxL + GxS + LxS+GxLxS + E, 
where Y is the observed value, fi is the overafi mean, 
G is tfie effect of tfie presence or absence of the P- 
element, 5 is the effect of sex, L is the random effect 
of the DGRP line versus CSB genetic backgrounds, 
GxL, Gx S, Lx S and GxLx S are the interaction 
terms, and e is the environmental variance between 
replicates. A significant interaction term (GxL and/ 
or Gx Lx S) indicates epistasis. To assess variation in 
epistatic effects among DGRP fines, we performed 
simfiar mixed model ANOVAs across afl genotypes, 
treating the DGRP genotypes and interactions with 
them as random effects. Finally, to determine the 
significance of epistatic interactions among different 
wild-derived genetic backgrounds, we estimated indi- 
vidual epistatic effects for each background and tested 
for significance using ANOVA. 
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3. Results 

(i) Effects o/Crc and Sema-5c mutants on olfactory 
behaviour 

To assess the effects of naturally segregating epistatic 
modifiers on P-element insertional mutations that 
affect olfactory behaviour, we selected P-element in- 
sertions in the Sema-Sc and Crc genes that have large 
effects on olfactory avoidance behaviour towards 
benzaldehyde (Sambandan et al., 2006; RoUmann 
et al., 2007). We verified the previously reported effects 
on responsiveness to benzaldehyde using our modified 
behavioural assay. To analyse the data, we used a 
two-way ANOVA model, Y=fi + L + S + LxS + E, 
where ft is the overall mean, L is the fixed effect of line, 
5" is the fixed effect of sex, L x 5 is the line x sex in- 
teraction term and E is the environmental variance. 

As we observed a significant line effect (P< 0-0001), 
but no significant sex (_P = 0-99) or line x sex effect 
(P=0-90), measurements of sexes separately were 
pooled for analyses. The RI at 0-3% (v/v) benzal- 
dehyde for the CSB control was 0-98 ±0-01 (« = 3 re- 
plicates/sex/genotype, 50 individuals per replicate), 
showing strong avoidance behaviour. RIs for the Sema- 
Sc and Crc mutants were 0-68 ±0-03 and 0-56 + 0 04, 
respectively (« = 3 replicates/sex/genotype, 50 in- 
dividuals per replicate for each mutant), significantly 
lower than the CSB control (/'<0-0001 ; Fig. 2a). 

(ii) Effects of Crc mutations on sleep 

Like Sema-5c, the Crc locus is a hotspot for /"-element 
insertions (Spradhng et al., 2011). Previously a P- 
element insertion allele at Crc (Crc^*^"^^**) was found 
to affect sleep phenotypes (Harbison & Sehgal, 2008), 
but this insertion was at a different site than the 
f-element insertion allele (Crc^^""^^) previously im- 
plicated in startle behaviour (Yamamoto et al., 2009). 
Although both insertions are in the first exon, it is 
possible that different insertion locations may have dis- 
tinct phenotypic effects (Rollmann et al., 2006, 2008). 
Therefore, we assessed the effects of the Qt^*^""^^ al- 
lele on day and night sleep and waking activity. There 
were significant differences between Crc^"^"™* and the 
co-isogenic control for night sleep in both sexes 
(/"< 0-0001), for day sleep for males (P< 0-0001), and 
for waking activity in females (P = 0-0024) (Fig. 2b-d). 
There was a significant sex effect for day sleep and 
waking activity (_P< 0-0001), and a significant sex x 
fine interaction for day sleep (/"< 0-0001, the mutation 
increases day sleep only in males). 

(iii) Epistasis between Crc attd Sema-5c mutations 
and wild-derived DGRP backgrounds 

We used inbred DGRP lines in which CJ has been 
replaced by a P-element free CSB wild-type Ci or a 
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Fig. 2. Effects of Crc and Sema-5c mutations on olfactory 
behaviour and sleep phenotypes compared with the 
co-isogenic CSB control, {a) Olfactory behaviour. Bars 
represent mean response indices for pooled sexes; error 
bars are standard errors of the mean, (b) Night sleep. 
(c) Day sleep, {d) Waking activity. Bars represent mean 
day and night sleep and waking activity for males and 
females, separately, for the CSB control (open bars) or the 
Crc mutant (black bars); error bars are standard errors of 
the mean. 
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Fig. 3. Observed (closed bars) and expected (open bars) mean response indices for olfactory behaviour of (a) 13 DGRP C3 
substitution lines with a P-element insertion at Sema-Sc and {b) 14 DGRP C3 substitution lines with a P-element insertion 
at Crc. The error bars indicate standard errors of the mean for pooled sexes, ns, not significant, *P<0-05, **P<0-01, 
***P<0 001. 



co-isogenic C3 carrying a /"-element insertion in Crc 
or Setna-Sc (Fig. 1 ; Yamamoto et at., 2009) to assess 
the effects of naturally segregating epistatic modifiers 
of the mutations on olfactory behaviour and sleep 
phenotypes. We measured olfactory behaviour for 13 
DGRP lines in which C3 was replaced by either an 
isogenic CSB C3 or a co-isogenic CSB chromosome 
with a Sema-5c mutation. We also assessed olfactory 
behaviour and sleep phenotypes for 14 pairs of DGRP 
fines with wild-type CSB and co-isogenic Crc C3. Since 
there were no significant effects of sex or sex x line 
interaction in the analyses of olfactory behaviour 
among the chromosome substitution fines, whereas 
these terms were significant for night sleep, day sleep 
and waking activity, we report the results for olfac- 
tory behaviour pooled across sexes, and the sleep and 
activity data for males and females separately. 

In the absence of epistasis, the expected phenotype 
of a DGRP line bearing a Senm-Sc or Crc mutation is 



the difference between the effect (2a) of the mutation 
in the CSB background and the observed mean pheno- 
type of the DGRP fine with the CSB C3. Epistasis is 
implicated by a significant difference between this ex- 
pected value and the observed mean phenotype of the 
DGRP line with a mutant C3. The significance of the 
estimate of the epistatic effect is given by the P-value 
of the genotype by fine interaction in an ANOVA 
comparing the effect of the mutation in CSB and the 
DGRP fine. Epistatic interactions that amplify the 
effect of the mutation are considered enhancer effects, 
whereas those that counteract the effect of the mu- 
tation are defined as suppresser effects. 

We found significant epistasis for olfactory behav- 
iour in all but one instance {Sema-Sc in RAL_437) 
(Fig. 3, Table 1). In all cases where significant 
epistatic interactions were observed for olfactory 
behaviour, the epistatic effects were negative; i.e. the 
observed responses of the substitution fines to 
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Table 1. Epistatic interactions for olfactory behaviour 
in DGRP chromosome substitution lines with Sema-5c 
or Crc mutations 



DGRP line 


Sema-5c 


Crc 


RAL 208 


-0-65*** 


-0-59*** 


RAL 303 


-0-33*** 


-0-26*** 


RAL 335 


-0-26*** 


-019** 


RAL 357 


nd 


Q.2Q*** 


RAL 358 


-0-33*** 


— 0-41*** 


RAL 360 


— 0-45*** 


nd 


RAL 362 


-0-15* 


0.^9*** 


RAL 365 


-0-46*** 


— 0-48*** 


RAL 375 


-0-56*** 


nd 


RAL 391 






RAL 399 


nd 


-0-34*** 


RAL 437 


-0-04 ns 


nd 


RAL 517 


nd 


-0-23** 


RAL 714 


nd 


— 0-42*** 


RAL 732 


Q.28*** 


-0-31*** 


RAL 786 


-0-36*** 


-0-19* 


RAL 852 


-0-31** 


-0-21** 



The values indicate estimated epistatic effects for olfactory 
RI of individual chromosome substitution lines with Sema- 
5c or Crc mutations. ***P<0-0001; **0-0001 <P<0-01 ; 
*0-01 <P<0-05; ns, P>0-05; nd, not determined. 

benzaldehyde were greater than predicted based on 
the estimate of 2a in the CSB background (Table 1, 
Fig. 3). Since the effect of the mutations is to reduce 
the response to benzaldehyde in the CSB background, 
the negative difference between observed and expec- 
ted olfactory behaviour in the DGRP hues indicates 
suppression of the mutant effect in wild-type back- 
grounds. We also found substantial and sex-specific 
epistasis between DGRP Knes and the Crc mutation 
for night sleep, day sleep and waking activity (Fig. 4, 
Table 2). For night sleep, epistatic interactions were 
mostly suppressing, as for olfactory behaviour, with 
few exceptions (e.g. RAL_358 and RAL_852 for 
females and RAL_365 for males; Fig. A a and h). The 
Crc mutation increases day sleep in males (Fig. 2b)\ 
thus, suppressing epistasis would counteract the Crc 
mutation by reducing day sleep duration. Interest- 
ingly, three genetic backgrounds showed epistatic 
interactions for day sleep for females, two of which 
were enhancer effects (Fig. 4c), indicating that muta- 
tions with no effects on a phenotype in one genetic 
background can have significant effects in other 
backgrounds (i.e. the effect of the mutation was sup- 
pressed in the CSB background). There were exten- 
sive epistatic effects for male day sleep (Fig. 4 J). 
These effects were exclusively suppresser effects; that 
is, epistasis caused day sleep duration to be dimin- 
ished in mutants that gave rise to prolonged day sleep. 
Few epistatic effects were observed for waking ac- 
tivity, with suppresser effects for both sexes in the 
RAL_365 background and enhancer effects for 



females in RAS_391 and males in RAL_517 (Fig. Ae 
and/). 

We assessed whether there was significant variation 
in epistasis among the wild-type and mutant CJ sub- 
stitution fines for each trait, as indicated by a signifi- 
cant genotype (wild-type versus mutant) by fine 
(DGRP fine) interaction in the ANOVA. This term 
was significant for all traits (Tables 3 and 4). Thus, 
there is variation in the extent to which natural var- 
iants modify mutational effects. 

(iv) Pleiotropic epistatic effects 

In addition to their effects on olfactory behaviour 
(Sambandan et al., 2006; Rollmann et al., 2007), the 
Sema-Sc and Crc mutations also show reduced startle 
behaviour (Yamamoto et al., 2008) and a mutation at 
Crc has been associated with reduced night and day 
sleep, and increased waking activity (Harbison & 
Sehgal, 2008). To assess whether the same epistatic 
modifiers affect the effects of the Sema-5c and Crc 
mutations on multiple traits, we first asked whether 
there was a correlation between the estimates of epi- 
static effects for olfactory behaviour and those of 
startle-induced locomotion, measured previously on 
the same lines (Yamamoto et al., 2009). We did not 
observe a significant correlation for either Setna-Sc 
(Fig. 5a) or Crc (Fig. 5b). Similarly, epistasis of 
olfactory behaviour was not significantly correlated 
with epistasis for day sleep and night sleep for Crc, 
and the correlation with waking activity in males was 
only nominally significant (^ = 0-04; Fig. 6). Epistasis 
of day time sleep, night time sleep or waking activity 
was also not correlated with epistasis of startle 
behaviour (Supplementary Fig. SI avafiable at http:// 
journals.cambridge.org/grh). These results show that 
different naturally segregating epistatic modifiers 
modulate different phenotypes affected by mutations 
at the same pleiotropic gene. 

4. Discussion 

Previously, olfactory behaviour in D. melanogaster 
has been used as a model trait to dissect the genetic 
architecture of behaviour (Anholt, 2010) and dynamic 
epistatic networks of pleiotropic genes have been im- 
plicated as a major feature of the genetic ensembles 
that underfie the manifestation of this behavioural 
phenotype (Fedorowicz et al., 1998; Sambandan 
et al., 2006). D. melanogaster can also serve as a gen- 
etic model to study sleep (Hendricks et al., 2000; 
Shaw et al., 2000). While epistasis can be hypothe- 
sized from co-regulated gene expression networks 
(Harbison et al., 2009), no previous study has quan- 
tified the impact of epistasis on sleep in flies. In ad- 
dition to modulation of behavioural phenotypes, 
suppressing epistasis may explain the paradox 
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Fig. 4. Observed (closed bars) and expected (open bars) sleep phenotypes in DGRP C3 substitution lines with a P-element 
insertion at Crc. (a) Night sleep in females, (h) Night sleep in males, (c) Day sleep in females, (d) Day sleep in males, 
(e) Waking activity in females. ( /') Waking activity in males. The error bars indicate standard errors of the mean for sexes 
separately, ns, not significant, *P<0 05, **P<0 01, ***P<0-001. 



between developmental robustness in the face of gen- 
etic variation, as illustrated by the elfects of genetic 
background modifiers on mutations in Sevenless and 
Drosophila Epidermal Growth Factor Receptor that 
affect development of photoreceptors (Polaczyk et al., 
1998). 

The recent generation of a panel of diverse homo- 
zygous wild-derived chromosome substitution fines 
that carry the same homozygous CSB C3 with or 



without a P-element insertion (Yamamoto et al., 
2009) enables analyses of the effects of naturally seg- 
regating epistatic modifiers. We used chromosome 
substitution lines with Sema-Sc and Crc mutations 
to analyse epistatic modulation of mutations that 
affect olfactory behaviour, sleep and waking activity 
in Drosophila in wild-derived genetic backgrounds. 
Sema-5c has been implicated in early development 
(Khare et al., 2000) and Crc, a calcium-binding 
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Table 2. Epistatic interactions for sleep phenotypes in DGRP chromosome 
substitution lines with a Crc mutation 



Waking activity 

DGRP line Sex Night sleep (min) Day sleep (min) (counts/min) 



RAL_208 


f 


-169-62 


* * * 


-24-46 


ns 


-0-23 


ns 




m 


-100-31 


* * * 


117-49 


* 


-0-33 


ns 


RAL_303 


f 


-19-19 


ns 


52-43 


ns 


-0-17 


ns 




m 


-37-75 


ns 


112-24 


* 


-0-18 


ns 


RAL_357 


f 


-96-62 


** 


-35-94 


ns 


0-18 


ns 




m 


-49-24 


ns 


224-17 


* 


0-32 


ns 


RAL_358 


f 


97-85 


* 


-42-37 


ns 


-0-21 


ns 




m 


108-11 


* * * 


47-16 


ns 


0-37 


ns 


RAL_362 


f 


-120-66 


* * * 


-66-80 


ns 


0-16 


ns 




m 


-59-56 


** 


122-30 


* ** 


-0-33 


ns 


RAL_365 


f 


73-22 


ns 


-199-98 


* ** 


-0-45 


** 




m 


162-15 


*** 


76-78 


* 


-0-72 


** 


RAL_391 


f 


-153-71 


*** 


125-73 


** 


0-46 


** 




m 


-63-63 


** 


233-48 


* 


0-21 


ns 


DAT "500 


f 
I 


— Z J(5 


ns 


A 1 .Q'^ 

— 01 yj 


ns 


n. 1 7 
— U 1 / 


ns 




m 


9-90 


ns 


17-48 


ns 


0-33 


ns 


RAL_517 


f 


-154-34 


* * * 


-87-52 


* 


-0-07 


ns 




m 


-115-57 


* * * 


91-86 


** 


0-50 


* 


RAL_714 


f 


-40-84 


ns 


-11-99 


ns 


-0-13 


ns 




m 


-24-10 


ns 


71-39 


* 


-0-28 


ns 


RAL_732 


f 


47-79 


ns 


13-18 


ns 


-0-01 


ns 




m 


51-55 


ns 


123-30 


* 


-0-11 


ns 


RAL_786 


f 


-173-88 


* * * 


-48-61 




-0-02 


ns 




m 


-132-34 


* * * 


116-61 


* 


0-08 


ns 


RAL_852 


f 


86-59 


** 


-91-57 


* 


-0-21 


ns 




m 


-55-06 


* 


14-67 


ns 


0-04 


ns 



The values indicate estimated epistatic elfects of individual chromosome substitution 
lines with a Crc mutation, m, males; f, females; *** P<0-0001 ; **0-0001 <P<0-01 ; 
*0 01 <P<0 05; ns, P>0-05. 

Table 3. ANOVAs of olfactory behaviour among DGRP lines with CSB and Sema-5c or Crc mutant third 
chromosomes 



Mutation 


Source of variation 


df 


SS 


MS 


F 


P 




Sema-Sc 


Genotype (G) 


1 


0-330 


0-330 


23-97 


< 0-0001 


* 




Line (L) 


12 


1-740 


0-102 


7-42 


< 0-0001 


** * 




Sex (S) 


1 


0-019 


0-019 


1-39 


0-24 


ns 




SxL 


12 


0-127 


0-007 


0-54 


0-9252 


ns 




GxL 


12 


0-929 


0-077 


5-62 


< 0-0001 


** * 




SxG 


1 


0005 


0-005 


0-36 


0-5495 


ns 




SxGxL 


12 


0-084 


0-007 


0-51 


0-9059 


ns 




Error 


104 


1-514 


0-014 








Crc 


Genotype (G) 


1 


0-078 


0-078 


6-37 


0-0129 


* 




Line (L) 


13 


2-468 


0-145 


11-92 


< 0-0001 


** * 




Sex (S) 


1 


0 000 


0-000 


001 


0-9206 


ns 




SxL 


13 


0-175 


0-010 


0-85 


0-6371 


ns 




GxL 


13 


0-599 


0-046 


3-78 


< 0-0001 


** * 




SxG 


1 


0-006 


0-006 


0-52 


0-4729 


ns 




SxGxL 


13 


0173 


0-013 


1-09 


0-3705 


ns 




Error 


112 


1-362 


0-012 









df, degrees of freedom; SS: sums of squares (type 111); MS, mean squares; ***/"< 0-0001 ; *0-01 <P<0-05; ns, not significant. 



chaperone, is involved in intracellular protein trans- 
port, exocytosis and development of the nervous 
system in Drosophila (Prokopenko et al., 2000). 
Mutations in Sema-5c reduce olfactory avoidance 



behaviour (Sambandan et al., 2006; Rollmann et al., 

2007) and startle behaviour (Yamamoto et al., 

2008) . Mutations in Crc result not only in aberrant 
chemosensory responses (Stoltzfus et al., 2003; 



Extensive epistasis for olfactory behaviour 17 



Table 4. ANOVAs of sleep phenotypes and waking activity among DGRP lines with CSB and Crc mutant third 
chromosomes 



Trait 


Source of variation 


df 


SS 


MS 


F 


P 




Night sleep 


Genotype (G) 


1 


828 359 


828 359 


7-40 


<0-0186 






Line (L) 


12 


5 493 950 


457 829 


2-10 


0-0654 






Sex (S) 


1 


73 388 


73 388 


0-62 


0-4462 


ns 




SxL 


12 


1426 496 


118875 


8-85 


< 0-0003 






GxL 


12 


1 349 197 


112433 


8-38 


< 0-0004 






SxG 


1 


6658 


6658 


0-50 


0-4942 


ns 




SxGxL 


12 


161 098 


13 425 


1-77 


0-0496 


ns 




Error 


TAT 
/U / 


5 ibo UUy 


/ jyz-oo 








Day sleep 


Genotype (G) 


1 


284 741 


284 741 


4-95 


< 0-0460 






Line (L) 


12 


5 543 930 


461 994 


4-69 


< 0-0002 






Sex (S) 


1 


4724 142 


4 724 142 


83-57 


< 0-0001 






SxL 


12 


681 402 


56 783 


3-53 


<0-0190 






GxL 


12 


693 050 


57 754 


3-59 


<0-0179 






SxG 


1 


1819 


1819 


0-11 


0-7423 


ns 




SxGxL 


12 


193 268 


16 106 


1-92 


0-0295 






Error 


707 


5 939 596 


8401-13 








Waking activity 


Genotype (G) 


1 


5-57 


5-57 


6-67 


0-0240 


* 




Line i^j^) 


1 0 
IZ 




I 50 


1 -AA 


0-2656 


ns 




Sex (S) 


1 


58-42 


58-42 


67-90 


< 0-0001 


** * 




SxL 


12 


10-37 


0-86 


2-11 


0-1051 


ns 




GxL 


12 


10-05 


0-84 


2-05 


0-1136 


ns 




SxG 


1 


0-00 


0-00 


0-00 


0-9660 


ns 




SxGxL 


12 


4-90 


0-41 


2-74 


0-0012 


** 




Error 


707 


105-57 


0-15 








df, degrees of 


freedom; SS, sums of 


squares (type 


III); MS, 


mean squares; 


***P<0-0001; 


**0-0001 <P<0-01; 



*0-01 <P<0-05; ns, not significant. 

Sambandan et al., 2006) and reduced startle behav- 
iotjr (Yamamoto et al., 2008), but also reduce day and 
night sleep duration and increase waking activity 
(Harbison & Sehgal, 2008). We confirmed the effects 
of these P-element insertions in the CSB background 
on olfactory behaviour using a recently developed 
modified high throughput olfactory behavioural assay 
(Swarup et al., 2011; Fig. 2a) and confirmed the ef- 
fects on sleep, using the same /"-element insertion fine 
in Crc previously implicated in startle behaviour 
(Yamamoto et al., 2009) and olfaction (Sambandan 
etai, 2006) (Fig. 2b). 

We found that mutational effects were generally 
reduced in the chromosome substitution lines com- 
pared with the original effect observed in the CSB 
background. The presence of variation in epistatic 
effects for each phenotype for each P-element inser- 
tion indicates that different wild-derived genetic 
backgrounds harbour different segregating epistatic 
modifiers that alter the effect of the P-element 
mutation. Although phenotypic measurements of a 
larger number of chromosome substitution fines 
might reveal correlations in epistatic measures among 
olfactory behaviour, startle behaviour and sleep, the 
lack of correlation of epistatic effects across these 
phenotypes among the 27 lines that were available for 
our study (Figs 5 and 6) suggests that different 
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_ 25 

1 ^° 
"i" 15 
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(b) Crc 
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0 0-2 0-4 0-6 ( 

Olfactoty behaviour (/□!,) 

Fig. 5. Relationship between the estimates of epistatic 
interactions for olfactory behaviour {Io\() and startle 
induced locomotion (/startle) in DGRP Ci substitution 
lines, (a) Seina-Sc: r^ = 0, />0-05. (h) Crc: /•2 = 0-186, 
P>0-05. 
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Fig. 6. Relationship between the estimates of epistatic interactions for olfactory behaviour (/oif) and sleep phenotypes 
(/sleep) in DGRP C3 substitution lines with a P-element insertion at Crc (a) Night sleep. Females: ;'^ = 0 001, P>0-05. 
Males: ;-' = 0-137, P>0-05. (b) Day sleep. Females: r^ = 0, P>0 05. Males: = 0-008, P>0 05. (c) Waking activity. 
Females: ;-' = 0-086, P>0-05. Males: r' = 0-319, P = 0-044. 




Fig. 7. Epistasis and pleiotropy. The diagram illustrates a 
focal P-element-tagged gene (red circle) that forms part of 
three genetic networks affecting different phenotypes, 
indicated by green, blue and orange colours, respectively. 
Gene ensembles that generate phenotype-specific epistatic 
interactions with the focal gene, indicated by the dotted 
arrows, are shown in corresponding muted colours. 



epistatic modifiers are likely to interact with the same 
pleiotropic gene to modulate different phenotypes 
(Fig. 7). This complex genetic architecture is in Hne 
with previous conclusions that the manifestation of 
complex behavioural phenotypes can be altered by 
ensembles of epistatic genes (Sambandan et al., 2006; 
Anholt, 2010; Zwarts et al., 2011). Independent 
segregation of components of these ensembles in a 
natural population will result in variation in epistatic 
effects and these effects may express themselves 
differently for different pleiotropic phenotypes as- 
sociated with the same causal variant. 

In conclusion, we have shown that epistasis appears 
to be a pervasive general feature of natural popu- 
lations and our results suggest that epistatic inter- 
actions may protect against adverse effects of new 
mutations. Furthermore, different epistatic inter- 
actions modulate different phenotypes affected by 
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mutations at the same pleiotropic gene. The 
prevalence of epistasis in the genetic architecture of 
complex traits is relevant to the design and interpret- 
ation of genetic studies in human populations. 
Widespread suppressing epistasis may account for the 
' missing heritability ' for human traits, such as height 
(Manoho et ai, 2009). Our study underscores the 
importance of D. melanogaster as a model system for 
the analysis of quantitative traits, as a similar detailed 
analysis of epistasis under conditions in which we can 
introduce a mutation in a range of tightly controlled 
genetic backgrounds would not be possible in human 
populations. Substitution of chromosomes with P- 
element insertions in DGRP backgrounds wiU enable 
future mapping of epistatic modifiers and, ultimately, 
genome-wide characterization of epistatic interactions 
between defined alleles and transposon-tagged sites 
that affect organismal phenotypes. 

This work was supported by grants from the National 
Institutes of Health (GM45146, GM59469) to TFCM and 
RRHA. 
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